— . . A  .  .  . .  . . .  — .  - 

^  REPORT  DOCUMENTATION  PAGE  AFRL-SR-AR-TR-06-0137 

Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  inst  3  data  n< 

and  completing  and  reviewing  this  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  o  D  Depar 

Defense,  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports  (0704-0188),  1215  Jefferson  Davis  Hie  aware  t 

notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  penalty  for  failing  to  comply  with  a  collection  of  information  i;  )  NOT  F 

YOUR  FORM  TO  THE  ABOVE  ADDRESS. _  _ _  _  _ 


1.  REPORT  DATE  (DD-MM-YYYY) 
10-March  2006 

2.  REPORT  TYPE 

Final  Performance  Report 

3.  DATES  COVERED  (From  -  To) 

1  Sep  2002  -  31  Oct  2005 

4.  TITLE  AND  SUBTITLE 

Mnlt-i.-PhagP.JJLeMq  nniming 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

F49620- 02 -1-0431 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

Suresh  Menon 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

8.  PERFORMING  ORGANIZATION  REPORT 

NUMBER 

1 

Georgia  Tech  Research  Corporation 

Georgia  Institute  of  Technology 

Atlanta,  GA  30332-0420 

9.  SPONSORING  /  MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES)  |  10.  SPONSOR/IV 

The  Air  Force  Office  of  Scientific  Research 

875  N.  Randolph  Street,  Suite  325,  Room  3112 

Arlington,  VA  22203-1768 

lONITOR’S  ACRONYM(S) 

AFOSR/NM 

11.  SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 

12.  Distribution  /  availability  statement 

Distribution  A:  Approved  for  Public  Release;  Distribution  Unlimited 

13.  SUPPLEMENTARY  NOTES 

14.  ABSTRACT 

One  Agent  Defeat  scenario  of  primary  interest  to  the  US  Air  Force  involves;  the  explosion  of  a  warhead  in  ground-fixed  target 
structures  and  storage  facilities.  Since,  improper  deployment  can  have  potentially  hazardous  consequences,  numerical  studies  can 
understand  the  dynamics  of  this  process.  Simulation  of  this  process  requires  a  multi  physics  code  capable  of  treating  both  shocks  £ 
turbulence  in  complex  geometries  as  well  as  single  and  multiphase  chemistry.  Smooth  flow  solvers  can  simulate  turbulent  flows,  l 
These  codes  usually  have  significant  difficulty  in  treating  discontinuities  such  as  shock  waves  without  introducing  unphysical 
oscillations  into  the  flow. 


20060601078 

15.  SUBJECT  TERMS 


16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION 

OF  ABSTRACT 

18. 

NUMBER 

19a.  NAME  OF  RESPONSIBLE  PERSi 

Suresh  Menon 

a.  REPORT 

Unclassified 

b.  ABSTRACT 

Unclassified 

C.  THIS  PAGE 

Unclassified 

. . .  . - 

UL 

19b.  TELEPHONE  NUMBER  (include  ar 

404-894-9126 

Standard  Form  298  (Rev.  8-98) 


Prescribed  by  ANSI  Std.  239.18 


SIMULATION  OF  SHOCK-SHEAR  INTERACTIONS  IN  COMPLEX  DOMAINS 

j  F49620-02-1-0431 

Suresh  Menon  and  Bruce  Fryxell 
Computational  Combustion  Laboratory 
School  of  Aerospace  Engineering 
Georgia  Institute  of  Technology 


Abstract 

One  Agent  Defeat  scenario  of  primary  interest  to  the  US  Air  Force  involves  the  explosion 
of  a  warhead  in  ground-fixed  target  structures  and  storage  facilities.  Since,  improper 
deployment  can  have  potentially  hazardous  consequences,  numerical  studies  can  help 
understand  the  dynamics  of  this  process.  Simulation  of  this  process  requires  a  multi¬ 
physics  code  capable  of  treating  both  shocks  and  turbulence  in  complex  geometries  as 
well  as  single  and  multiphase  chemisty.  Smooth  flow  solvers  can  simulate  turbulent 
flows,  but  these  ;odes  usually  have  significant  difficulty  in  treating  discontinuities  such 
as  shock  waves!  without  introducing  unphysical  oscillations  into  the  flow.  Shock 
capturing  methods  do  a  good  job  at  treating  shocks  and  contact  discontinuities,  but  the 
dissipation  in  thdse  methods  can  provide  an  incorrect  treatment  of  the  turbulent  features 
in  the  flow.  In  this  project,  we  have  developed  a  hybrid  code,  which  combines  a  high- 
order  smooth-fldw  solver  with  a  high-resolution  shock-capturing  method.  The  shock- 
capturing  method,  is  used  only  near  discontinuities,  while  the  smooth-flow  solver  is  used 
to  treat  the  remajnder  of  the  flow.  The  hybrid  method  is  capable  of  producing  solutions 
that  are  more  accurate  than  could  be  obtained  with  either  method  by  itself. 

The  baseline  code  is  LESLTE3D1’2’3’4,  a  well-established  DNS/LES  solver  developed  for 
scalar  mixing  arid  combustion  applications.  The  code  contains  a  smooth-flow  Navier- 
Stokes  solver  that  is  second-order  accurate  in  time  and  fourth  order  in  space.  It  employs 
a  localized  dynamics  model  for  the  subgrid  closure  of  the  unresolved  momentum  and 
energy  fluxes  thht  has  been  shown  to  be  highly  robust  and  stable  in  all  types  of  flows. 
There  are  no  adjustable  constants  or  coefficients  in  this  model,  since  the  coefficients 
adapt  automatically  to  the  problem  of  interest.  The  code  has  the  ability  to  deal  with 
multiphase  (gas-liquid-solid)  flows  using  a  combined  Eulerian-Lagrangian  approach,  in 
which  the  gas  phase  is  treated  on  a  fixed  Eulerian  grid,  while  the  liquid  and  solid  particles 
are  tracked  in  a  grid-free  Lagrangian  manner.  The  code  has  been  extensively  verified  and 

validated  against'a  wide  variety  of  experimental  and  canonical  flows. 

| 

Although  LESL1E3D  works  very  well  for  smooth  flows,  it  cannot  accurately  handle 
strong  shocks  arid  contact  discontinuities,  which  produce  unphysical  oscillations  in  the 
flow.  A  major  accomplishment  under  this  grant  has  been  to  incorporate  a  high-resolution 
shock-capturing  [method  into  the  code.  The  scheme  we  use  for  shock  capturing  is  the 
P i cce wi se-Parabolic  Method  (PPM)  of  Colelia  and  Woodward5.  This  method  has  been 
widely  used,  particularly  for  astrophysics  calculations,  where  strong  shocks  appear  on  a 
regular  basis,  the  implementation  used  here  is  based  on  the  PROMETHEUS  code6, 
which  is  also  the1  basis  of  the  hydrodynamic  solver  in  the  FLASH  code7,  developed  at  the 


University  of  Chicago  as  part  of  the  DOE  ASCI  Alliance  Program.  The  code  has  been 
used  by  many  groups  throughout  the  world  and  has  been  extensively  validated8'9. 

PPM  is  a  high-order  extension  of  Godunov's  method10.  The  algorithm  was  designed 
specifically  to  produce  accurate  solutions  of  flows  with  shaip  discontinuities  through  the 
use  of  nonlinearity  and  ''smart"  dissipation  algorithms.  It  also  performs  very  well  in 
smooth  flows.  Godunov’s  approach  is  to  calculate  the  fluxes  in  the  gas  dynamic  equations 
by  solving  Rieminn's  problem  at  each  zone  interface.  This  produces  an  upwind-centered 
flux  and  also  injroduces  nonlinearity  into  the  difference  equations.  This  nonlinearity 
permits  the  method  to  calculate  narrow  discontinuities  without  generating  spurious 
oscillations.  Godunov's  method  assumes  a  pieccwise-constant  profile  of  each  variable. 

This  limits  the  accuracy  of  the  scheme  to  first  order. 
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A  formalism  for!  extending  Godunov's  method  to  higher  order  was  developed  with  the 
MUSCL  scheme!1,  in  which  second-order  accuracy  in  both  space  and  time  is  achieved  by 
representing  the;  flow  variables  as  piecewise-linear  instead  of  piecewise-eonstant 
functions.  Another  major  advance  of  MUSCL  was  the  incorporation  of  monotonicity 
constraints  instead  of  artificial  viscosity  to  eliminate  oscillations  in  the  flow.  PPM  takes 
the  next  logical  Step  by  representing  the  flow'  variables  as  piecewise-parabolic  functions. 
This  corresponds!  to  switching  to  Simpson's  rule  for  numerical  integration. 

The  problem  of  aj  point  explosion  in  a  confined  domain  was  performed  as  an  initial  test  of 
the  accuracy  and;  robustness  of  the  shock-capturing  method  for  agent  defeat  simulations. 
The  calculation  (;ests  the  ability  of  the  code  to  capture  both  strong  and  weak  shocks 
accurately  and  its  ability  to  compute  complex  shock  reflections  and  interactions.  The 
initial  setup  consists  of  a  cubical  box  25  m  wide  containing  air  at  room  temperature.  An 
explosion  is  initiated  off-center  in  the  box  by  raising  the  pressure  at  a  single  grid  point  by 
a  factor  of  106.  Figure  1  shows  the  results  of  a  2D  600  x  600  simulation  at  four  different 
times.  The  boundaries  were  treated  as  slip  walls.  The  explosion  w?as  initiated  at  a  point 
1/3  of  the  way  across  the  box  from  the  left  and  bottom  walls.  As  the  simulation  evolves, 
the  reflections  of  the  blast  wave  off  the  U'alls  of  the  box  and  the  resulting  shock 
interactions  resu.t  in  a  very  complex  pattern  of  waves.  The  shock  strength  rapidly 
decays,  and  by  the  end  of  the  simulation,  the  maximum  overpressure  is  less  than  1% 


Figure  1 .  Time  sequence  of  the  pressure  distribution  resulting  from  a  point  explosion  in  a 
box.  The  simulation  was  performed  with  PPM  in  two  dimensions  using  a  300  x  300  grid. 


A  three-dimensional  simulation  of  the  problem,  performed  on  a  192  x  192  x  192  grid  is 
shown  in  Figure  12.  For  this  simulation,  a  small  hole  was  inserted  into  one  wall  to  allow 
some  of  the  explosion  energy  to  escape. 


Figure  2.  Time  sequence  of  the  pressure  distribution  resulting  from  a  point  explosion  in  a 
box.  The  simulation  was  performed  with  PPM  in  three  dimensions  using  a  1923  grid. 

In  this  study,  wji  have  developed  a  new  hybrid  method,  which  uses  PPM  only  near 
discontinuities  and  calculates  the  remainder  of  the  flow  using  the  smooth-flow  solver. 

We  define  a  smoothness  parameter  given  by 

! 

s  \0.i-2Q+Qm\ 

where  <2;  is  any  variable  of  interest.  This  is  similar  to  the  refinement  criterion  for 
adaptive  mesh  refinement  proposed  by  L6hner12.  Typically,  we  use  both  pressure  and 
density.  To  avoid  detecting  numerical  noise,  the  smoothness  parameter  is  set  to  0  if 
either  the  numerator  or  denominator  is  less  than  some  small  value,  typically  0.05.  For 
multidimensional  simulations,  we  take  the  maximum  of  the  smoothness  parameter  for 
each  coordinate  direction,  ignoring  cross  derivatives.  If  the  smoothness  parameter  is 
greater  than  0.35  at  a  grid  point,  we  use  the  PPM  fluxes  to  update  the  conserved 
quantities  at  that  point.  Where  the  smoothness  parameter  is  less  than  0.35.  the  fluxes  from 
the  smooth-flow  [solver  are  used.  In  order  to  minimize  any  numerical  noise  that  might  be 
generated  when  [switching  between  the  two  methods,  we  take  the  average  of  the  two 
fluxes  at  any  grid  point  bordering  a  point  in  which  the  PPM  fluxes  are  needed. 

Wc  illustrate  the  hybrid  method  by  simulating  Richtmyer-Meshkov  instability'.  The  left 
side  of  the  grid  contains  hydrogen,  while  the  right  side  of  the  grid  consists  of  helium.  A 
sinusoidal  interface  is  placed  between  the  two  fluids,  with  a  second  higher-order 
perturbation  with  1/5  the  wavelength  and  1/10  the  amplitude  superimposed  on  the  large- 
scale  sinusoidal  perturbation.  A  planar  shock  with  an  overpressure  of  300  is  place  just  to 
the  left  of  the  interface.  The  results  of  a  two  dimensional  simulation  on  a  1024  x  256  grid 
arc  shown  in  Figure  3  and  the  3D  simulation  using  a  480x120x120  grid  is  shown  in  Fig. 
4.  The  instability  shows  the  characteristic  bubble  and  spike  morphology.  The  shock  is 
located  just  to  the  right  of  the  left  boundary.  A  diamond-shaped  pattern  of  contact 
discontinuities,  which  arc  subject  to  Kelvin-Helmholtz  instabilities,  develops  between  the 
shock  and  the  interface,  and  a  series  of  Kelvin-Helmholtz  instabilities  forms  along  the 
interface.  The  tep  frame  in  Fig,  3  shows  the  magnitude  of  the  density  gradient  obtained 
with  PPM,  the  middle  frame  shows  the  results  of  the  hybrid  method,  and  the  bottom 
frame  shows  where  each  flux  was  used.  The  dark  gray  regions  were  calculated  with  the 
PPM  fluxes,  the  light  gray  regions  used  the  fluxes  from  the  smooth-flow  solver,  and  the 
average  of  the  two  fluxes  was  used  in  the  black  regions.  PPM  and  the  hybrid  method 
give  virtually  identical  results. 


The  hybrid  method  is  also  being  used  for  LES  of  turbulence.  The  simulation  shown  here 
is  of  a  planar  shock  propagating  through  isotropic  turbulence.  The  initial  conditions  were 
taken  from  a  DNS  simulation  of  Mahesh.  Lele  and  Moin13,  w'ho  used  a  6th  order  method 
on  a  231  x  81  x  81  grid.  Our  simulation  was  performed  with  the  hybrid  method  on  a  62  x 
32  x  32  grid.  Th;e  average  turbulent  kinetic  energy  as  a  function  of  distance  is  shown  in 
Figure  5.  The  plot  on  the  left  shows  the  results  obtained  using  the  smooth-flow  solver  by 
itself  and  using:  PPM  by  itself.  The  diamonds  represent  the  results  of  the  DNS 
simulation.  The; smooth-flow  solver  (solid  line)  produces  a  shock  that  is  too  wide  and 
overestimates  the  amount  of  post-shock  turbulence.  PPM  underestimates  the  peak  of  the 
post-shock  turbujence.  The  plot  on  the  right  shows  the  result  of  the  hybrid  method.  This 
result  matches  the  results  of  the  DNS  very  well. 
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The  final  simulation  presented  here  is  a  ID  multi-phase  detonation.  The  detonation  is 
ignited  in  a  mixture  of  llh  +  O2+  3 Ar,  originally  at  26  kPa.  The  particles  are  1 5  pm  A1 
particles  with  a  mass  loading  of  20  g/m3.  A  7-step,  7-species  H2-02  mechanism  is  used. 
A1  combustion  is  by  a  3-step  reaction  mechanism.  Pressure  traces  at  twm  points 
downstream  from  the  ignition  point  are  shown  in  Figure  6.  The  A1  particles  absorb  some 
of  the  reaction  energy,  heating  up  until  they  start  evaporating.  The  heat  exchange 
between  the  pariiclcs  and  the  gas  acts  as  an  energy  loss  to  the  gas,  thus  generating 
pressure  losses  ill  the  post-shock  region.  As  a  consequence,  the  propagation  velocity  of 
the  detonation  is  slower  when  the  particles  are  included. 


Figure  3.  Two-dimensional  simulation  of  a  R-M  instability  on  a  1024  x  256  grid. 


Figure  4.  Temperature  iso-surfaces  from  a  three-dimensional  simulation  of  a  Richtmyer- 
Meshkov  instability  computed  using  the  hybrid  method  on  a  480  x  120  x  120  grid. 


The  final  simulation  presented  here  is  of  a  one-dimensional  multi-phase  detonation.  The 
detonation  is  ignited  in  a  mixture  of  2H2  +  O2  +  3  Ar,  originally  at  26  kPa.  The  particles 
used  in  the  simulation  were  15  pm  A1  particles  with  a  mass  loading  of  20  g/m3.  The 
chemistry  of  thejgaseous  phase  is  treated  with  a  7-step,  7-species  mechanism.  Once  tire 
A1  particles  evaporate,  they  are  allowed  to  react  with  the  gaseous  phase  using  a  3-step 
reaction  mechanism.  Pressure  traces  at  two  points  downstream  from  the  ignition  point 
arc  shown  in  Figure  6  for  detonations  both  with  and  without  the  A1  particles.  The  Al 
particles  absorb  some  of  the  reaction  energy,  heating  up  until  they  start  evaporating.  The 
heat  exchange  between  the  particles  and  the  gas  acts  as  an  energy  loss  to  the  gas,  thus 
generating  pressure  loss  in  the  post-shock  region.  As  a  consequence,  the  propagation 
velocity  of  the  detonation  is  slower  when  the  particles  are  included. 


Figure  5.  Turbulent  kinetic  energy  obtained  from  an  LES  simulation  of  the  propagation 
of  a  planar  shock  through  isotropic  turbulence.  The  plot  on  the  left  shows  the  results 
obtained  with  tljc  smooth-flow  solver  (solid  line)  and  with  PPM  (dashed  line).  The 
diamonds  are  the  results  obtained  from  a  DNS  simulation. 
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time  from  ignition  (s) 


Figure  6.  Pressure  profiles  0.7  rn  and  1.4  m  downstream  from  the  ignition  point  for 
detonations  both  with  and  without  evaporating  and  reacting  Al  particles. 


Acknowledgmeift/Disclaimer 

This  work  was  sponsored  (in  part)  by  the  Air  Force  Office  of  Scientific  Research,  USAF, 
under  grant/contfact  number  F49620-  F49620-02-1-0431.  The  views  and  conclusions 
contained  herein!  are  those  of  the  authors  and  should  not  be  interpreted  as  necessarily 
representing  the  official  policies  or  endorsements,  either  expressed  or  implied,  of  the  Air 
Force  Office  of  Scientific  Research  or  the  U.S.  Government. 

! 

References 

1 .  S.  Menon,  Computational  and  Modeling  Constraints  for  Large-Eddy  Simulations  of 
Turbulent  Combustion,  International  Journal  of  Engine  Research,  2000, 1:  209-227. 

2.  W.-W.  Kim  ajid  S.  Menon.  Numerical  Modeling  of  Turbulent  Premixed  Flames  in  the 
Thin-Gas  Turbine  Combuster,  Combustion  Sci  and  Technology,  1999, 143: 25-62. 

3.  W.-W.  Kim  dnd  S.  Menon,  Numerical  Modeling  of  Fuel-Air  Mixing  in  a  Dry'  Low- 
Reaction-Zorjes  Regime,  Combustion  Science  and  Technology,  1999,  60:  113-150. 

1 

4.  C.  Stone  and  S.  Menon,  Parallel  Simulations  of  Swirling  Turbulent  Flames,  Journal  of 
Supercomputjng,  2003,  22:  7-28. 

5.  P.  Colella  arid  P.  Woodward,  The  Piecewisc-Parabolic  Method  (PPM)  for  Gas 
Dynamical  Simulations,  Journal  of  Computational  Physics,  1984,  54:  174-201. 

6.  B.  Fryxcll,  E.  Miiller,  and  D.  Arnett,  Hydrodynamics  and  Nuclear  Burning,  Max- 
Planck-Institiite  fur  Astrophysik  Report  449,  1989. 

7.  B.  Fryxell,  K,-  Olson,  P.  Ricker.  F.X.  Timmes,  M.  Zingalc,  D.Q.  Lamb,  P.  MacNeice, 
R.  Rosner,  J.W.  Truran,  and  H.  Tufo,  FLASH:  An  Adaptive  Mesh  Hydrodynamics 
Code  for  Modeling  Astrophysica!  Thermonuclear  Flashes,  Astrophysical  Journal 
(Supplement),  2000,  131:  273-334. 

8.  P.  Woodward  and  P.  Colella,  The  Numerical  Simulation  of  Two-Dimensional  Fluid 
Flow  with  Strong  Shocks,  Journal  of  Computational  Physics,  1984,  54:  1 15-173. 

9.  A.  Calder,  B.  Fryxell,  T.  Plewa,  R.  Rosner,  L.  Dursi,  V.  Weirs,  T.  Dupont,  H.  Robey, 
J.  Kane,  B.  Remington,  R.  Drake,  G.  Dimonte,  M.  Zingale,  F.  Timmes,  K.  Olson,  P. 
Ricker,  P.  MacNeice,  and  H.  Tufo,  On  Validating  an  Astrophysical  Simulation  Code, 
Astrophysica  Journal  (Supplement),  2002, 143:  201-229. 

10.  S.  Godunov.  A  Finite  Difference  Method  for  the  Computation  of  Discontinuous 
Solutions  of  the  Equations  of  Fiuid  Dynamics,  Mat,  Sb.,  1959,  47:  357-393. 

11.  B.  van  Leer,  Towards  the  Ultimate  Conservative  Difference  Scheme  V.  A  Second 
Order  Sequel  to  Godunov's  Method,  Journal  of  Comp.  Physics,  1979,  32:  101-136. 

12.  R.  Lohner,  Computational  Methods  in  Applied  Mech  Engineering,  1987,  61:  323. 


13.  K.  Mahesh.  js.  Lclc,  and  P.  Moin,  The  Influence  of  Entropy  Fluctuations  on  the 
Interaction  of  Turbulence  With  a  Shock  Wave,  Journal  of  Fluid  Mechanics,  1997, 

334:  225-241 

| 

j 

Personnel  Supported  During  Duration  of  Grant 

Jason  Hackl  Graduate  Student,  Georgia  Institute  of  Technology',  Atlanta 

Bruce  Fryxell  Senior  Research  Scientist,  Georgia  Institute  of  Technology, 

Atlanta 

Suresh  Menon  Professor,  Georgia  Institute  of  Technology',  Atlanta 

Publications 

B.  Fryxell,  S.  Menon,  F.  Genin,  and  J.  Hackl,  A  Hybrid  Method  for  Large  Eddy- 
Simulations  With  Strong  Shocks,  AIAA  Journal,  to  be  submitted, 

Fryxell,  B.  and  Menon,  S., "Hybrid  Simulations  of  Riehtmycr-Meshkov  Instability," 
AIAA-2005-03 14, 43rd  AIAA  Aerospace  Sciences  Meeting,  Reno,  NV,  January,  2005. 

Menon,  S.,  Fryxell.  B.  and  Genin,  F.,  "Parallel  Simulations  of  Detonation  in  Multi-Phase 
Turbulent  Reacti  ig  Mixtures,"  2005  DOD  HPC  Users  Group  Conf.,  June  28, 2005. 

Genin,  F.,  Fryxell,  F,  and  Menon,  S.,  "Simulation  of  Detonation  Propagation  in 
Turbulent  Gas-Solid  Reactive  Mixtures, "AIAA-2005-3967,  42nd  AIAA/ASME/ASEE 
Joint  Propulsion  Conference,  July  11-14,  2005. 

Genin,  F.,  Fryxell,  B.  and  Menon,  S.,  "Hybrid  Large-Eddy  Simulation  of  Detonation  in 
Reactive  Mixtures,"  to  be  presented  at  the  20th  International  Conference  on  Detonations, 
Explosions  and  Shock  Waves,  Montreal  CA,  August  1-4, 2004. 

Invited  Talks 

Bruce  Fryxell,  Calculation  of  Turbulent  Flows  with  Strong  Shocks,  Los  Alamos  National 
Laboratory,  May  25, 2005 

AFRL  Point  of  Contact 

Douglas  Nance,  AFRL/MNAC,  Eglin  AFB,  FL,  Phone  850-882-8302  x  3208.  Regular 
bi-monthly  meeting  at  Georgia  Institute  of  Technology. 

Transitions 

Tire  LESLIE3D  and  LESLIE-PPM  codes  have  been  delivered  to  AFRL/MNAC  for  their 
in-house  classified  studies  of  the  Agent  Defeat  problems.  Successful  application  of  these 
codes  to  practical  problems  has  been  demonstrated  at  Eglin  AFB. 

New  Discoveries 


